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Abstract: While conventional hydrodynamics incorporating dissipative effects is 
hard to derive from an action principle, it is nevertheless possible to construct classi- 
cal actions when the dissipative terms are switched off. In this note we undertake a 
systematic exploration of such constructions from an effective field theory approach 
and argue for the existence of non-trivial second order non-dissipative hydrodynamics 
involving pure energy- momentum transport. We find these fluids to be character- 
ized by five second-order transport coefficients based on the effective action (a three 
parameter family is Weyl invariant). On the other hand since all flows of such fluids 
are non-dissipative, they entail zero entropy production; one can therefore under- 
stand them using the entropy current formalism which has provided much insight 
into hydrodynamic transport. An analysis of the most general stress tensor with 
zero entropy production however turns out to give a seven parameter family of non- 
dissipative hydrodynamics (a four parameter sub-family being Weyl invariant). The 
non-dissipative fluids derived from the effective action approach are a special case of 
the fluid dynamics constrained by conservation of the entropy current. We specu- 
late on the reasons for the mismatch and potential limitations of the effective action 
approach. 
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1 Introduction 

Fluid dynamics is an effective description, valid over spatio-temporal length scales 
large compared to the mean free path of any interacting system. Traditionally the 
variables of fluid dynamics are the local fluid velocity, temperature (or energy density) 
and the local charge densities. The dynamical equations are simply the conservation 
equations for the stress tensor and other conserved currents. The basic phenomeno- 
logical input of fluid dynamics are the constitutive relations where the stress tensor 
and other conserved currents are expressed in terms of the fluid variables i.e., velocity, 
temperature and the charge densities. 

But these constitutive relations cannot be completely arbitrary. Because fluid 
dynamics is always a corse-grained description of some underlying microscopic quan- 
tum theory, it must satisfy all the consistency requirements arising from the quantum 
dynamics of the fundamental degrees of freedom (in an appropriate long-wavelength 
fluid limit). 

First of all the constitutive relations are constrained by the symmetries of the 
system. But more importantly we must use a local version of second law to further 
constrain the phenomenology of the stress tensor and currents. In particular, we 
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demand that the entropy should locally increase on every time dependent solution 
of the fluid equations i.e., there should exist at least one entropy current whose 
divergence is always non-negative for every fluid flow consistent with the equations 
of motion. It turns out that this particular condition constrains the possible form of 
the constitutive relations to a large extent, a fact that has been well-appreciated for 
a while, cf., [1]. 

However, in spite of decades of work on fluid dynamics, we do not yet know 
if the aforementioned requirements suffice to guarantee the full consistency of the 
hydrodynamical equations, when we further demand that such a fluid is an infra-red 
effective description of some microscopic quantum theory. We do not yet have a 
formulation which can capture all these microscopic constraints simultaneously (nor 
have we established what the full set of constraints are). What we do know is that it 
appears to be sensible to use a working definition of hydrodynamics in terms of the 
symmetries and a local version of the second law, cf., [2-7] for recent developments 
which explore this philosophy in detail. 

Generically fluid dynamics involves dissipation; to leading order in gradients this 
is reflected by the parabolic nature of the energy-momentum conservation V^T^ = 
arising from viscous terms. This has entailed that traditional analyses of fluid 
dynamics as an effective field theory involves working with the equations of motion. 
It is difficult to write some quantity like an 'effective action' which upon extremization 
produces these equations thereby ensuring all such consistencies at one go. 

It is however intriguing to ask if it is possible to consider purely non-dissipative 
fluids. These clearly have a chance of being derivable from an effective action, and 
provide an opportunity to explore the microscopic constraints on fluid dynamics. 
Clearly, ideal hydrodynamics (a perfect fluid) is an example of such non-dissipative 
dynamics. A-priori it is not clear if we can deform away from perfect fluidity by 
higher order terms; if this were not possible then we would be discussing a rather 
boring class of examples. Our program of using effective actions will be essentially 
re-deriving thermodynamics (which clearly can be done using an Euclidean path 
integral to compute the free energy). Note that since we are switching off dissipative 
terms, we are guaranteed to have an entropy current which is identically conserved. 

Fortunately, it turns out that the space of non-dissipative fluid dynamical theo- 
ries is more than ideal fluids (as we argue below). One can proceed either by starting 
with a dissipative fluid system, setting the dissipative transport terms to zero and 
checking that the resulting system satisfies known constraints of an autonomous hy- 
drodynamics, viz., symmetry principles and entropy conservation. However, we can 
also proceed by trying to construct and effective action, since any action formal- 
ism as argued earlier could capture the physics of non-dissipative fluids. The main 
question is what is the organizing principle for such an action. It turns out the 
one can construct an effective action for hydrodynamics in rf-dimensional spacetime 
in terms of d — 1 scalar fields (ft 1 , which are the (Lagrangian) labels for local fluid 
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element as a function of the fixed spacetime coordinates x M . Alternately one can 
view the system in an Eulerian language by prescribing the spatial position as a 
function of the fluid element labels and time. The effective action is constructed as a 
re-parameterization invariant action respecting volume preserving diffeomorphisms 
in the space. This approach has been used to describe ideal fluids over the years 
motivated by an attempt to write action principles; cf., [8-17]. 

The most modern analysis of such effective actions for hydrodynamics was car- 
ried out recently in [16]. In particular, they argued that such an effective field the- 
ory naturally implements entropy conservation as an identity (a consequence of the 
volume preserving diffeomorphisms) thus making them ideal for explorations of non- 
dissipative hydrodynamics. For both uncharged and charged fluids they have shown 
how their formulation produces ideal fluid dynamics and they discuss some higher 
order corrections to exemplify their approach. 1 We will be interested in a general 
analysis of such classes of actions; we will classify all possible higher order corrections 
(to second order in an appropriate derivative counting scheme) and understanding 
the resulting conditions on the transport coefficients. 

Before proceeding however, we should first issue a disclaimer: while we will 
analyze the constraints arising on hydrodynamics from the existence of a classical 
effective action, we will not be able to argue that the resulting classical theory itself 
arises from an underlying unitary quantum dynamics upon coarse-graining. For much 
of the text we will analyze the classical description of non-dissipative hydrodynamics, 
relegating to §5 some speculations on whether one can realize such from well behaved 
quantum dynamics. 

The first task we undertake is to examine the effective action built from the 
Lagrangian fluid element variables systematically in a gradient expansion. We need at 
every order in the expansion to classify the number of on-shell inequivalent operators 
built out of the variables <f) ! . As mentioned earlier [16] analyze the leading order piece 
in such a derivative expansion (we explain our derivative counting more carefully 
below). The number of allowed structures grows with increasing order, but it easy 
to enumerate the relevant operators at two orders beyond the leading one. Focussing 
specifically on uncharged fluids, we find that the first non-zero correction to the 
action arises only at second order, where 5 distinct on-shell inequivalent operators 
exist. As a result we can have 5 free parameters 2 (see Eq. (2.15)) entering in our 
effective action and therefore generates a stress tensor with 5 independent transport 
coefficients (see Eq. (2.17)). 

lr rhis effective action formalism has been shown to reproduce certain features of anomalous 
hydro dynamic transport [17] in two spacetime dimensions as well as describing features of parity 
violation in 3 spacetime dimensions [15]. We will comment on these constructions at appropriate 
stages in our discussion. 

2 In this note, by the word 'parameter', we always mean an arbitrary function of temperature or 
entropy density. Thus our parameters always refer to transport functions. 
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We then analyse how the existence of an entropy current with zero divergence 
constrains the possible form of the stress tensor by generalizing the analysis of [7]. 
This reveals that one has a second order stress tensor with 7 independent transport 
coefficients cf., Eqs. (3.1), (3.12) below. 

Therefore the two stress tensors, derived using these two methods almost agree 
except these two free transport coefficients. We find that the stress tensor determined 
from action is a special case of the stress tensor determined using the entropy current 
analysis, provided we do some identification of the parameters, cf., Eq. (4.3). 

This is what we should expect since the 'action formalism' always admits a 
conserved entropy current by construction. However it is interesting that from this 
action we do not get the most general stress tensor allowed by the existence of a 
conserved entropy current. This can be interpreted either as an indication that the 
'entropy current technique' is not sufficient or on the contrary as a hint that some 
further generalization of this 'action formalism' is possible. 

This note is organized as follows. In §2 we shall briefly discuss the effective action 
formulation of non-dissipative fluids following the treatment of [16] and extend it 
to second order to derive the stress tensor. In §3 we shall specialize the entropy 
current analysis of [7] to the case of zero divergence entropy current and derive the 
constraints on the transport coefficients. In order to compare the two computations in 
§4 we first briefly describe the frame invariant formulation which allows to compare 
fluid dynamical data derived in inequivalent frames. We then show that one can 
sensibly compare the two stress tensors (and entropy currents) derived in §2 and 
§3 respectively and present the identification of the parameters appearing in the 
preceding two sections. In section §5 we shall conclude with some implications of 
our discussion and some open questions. We collect some useful intermediate results 
and conventions in Appendix A and present a brief analysis of parity violating fluids 
in 3 spacetime dimensions in Appendix B. 

2 Effective actions for non-dissipative hydrodynamics 

As mentioned in the introduction, it should be possible to construct and effective 
action for non-dissipative hydrodynamics. Such a construction should satisfy two 
primary requirements: 

• The Euler-Lagrange equations arising from the effective action which lead to 
the equations of motion of the theory should have no more content than energy- 
momentum conservation. The latter arise as a consequence of diffeomorphism 
invariance of the action, since the stress-tensor can be obtained by varying the 
action with respect to the background metric. 

• Lack of dissipative effects means that in such a formalism one should be able 
to identify a conserved current, which we interpret as the entropy current Jg. 
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The requirement of equations of motion having no more dynamical content than 
stress tensor conservation is a strong one. Usually while it is true that equations of 
motion imply V^T^ = the converse often fails to hold. 

The key point is identifying the correct degrees of freedom which achieves this; 
such a construction has been known for a long time based on what we might call, 
Lagangian fluid variables. Intuitively we want some symmetry in the field space 
that mimics the diffeomorphism invariance of the background spacetime. This can 
trivially be achieved by demanding that that configuration space of our classical 
action be parameterized by canonical field variables which respect field redefinition 
invariance (which is the analog of diffeomorphism in field space). Then by passing to 
some gauge fixed version a la,, static gauge we can argue that field variations which 
lead to Euler-Lagrange equations can be conflated with background diffeomorphisms 
thereby ensuring that energy-momentum conservation being the dynamical equations 
of the theory; cf., [18] for observations relating to this point. 

The fields of interest are labels for individual fluid elements viewed as a function 
of background spacetime coordinates. Such formulations have been described for 
perfect fluids (both relativistic and non-relativistic) for a very long time as mentioned 
in the introduction. We will use the recent analysis of [16] who were able to use these 
variables to motivate an effective action for non-dissipative hydrodynamics. While 
[16] discuss both neutral and charged fluids we will focus exclusively on neutral fluids 
in our discussion below. 

2.1 The fundamental fields of hydrodynamics 

Let us quickly review the ingredients in the construction of [16] starting with un- 
charged fluids. As mentioned earlier we want to work with local fluid elements, and 
use fields (ft 1 describe the position of the local fluid element in space at an instant of 
time. We will work with (i-dimensional fluids, so i = 1, 2, • • • , d — 1. We further fix 
the geometry of the spacetime in which the fluid propagates and take coordinates x M 
to be an appropriate chart on this spacetime manifold, with metric g^ u . 

Since we are trying to tag local fluid elements, we expect that the description 
of the low energy effective dynamics enjoys translational and rotational invariance 
of these co-moving coordinates. This tells us that the effective action should be 
constructed out of the derivatives of the fields (j) 1 and be suitably rotationally invari- 
ant. In addition, a local version of Liouville theorem in this configuration space of 
the cf) 1 demands invariance under arbitrary reparameterizations of the (j) 1 , i.e., the 
Lagrangian should be invariant under 

¥ -»■ £ J (0) , Jacobian(£, <p) = 1 (2.1) 

with the condition of the Jacobian being forced on us by the fact the volume of 
configuration space be unchanged. Clearly, such a symmetry is generated by the 
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diffeomorphisms of the {(f) 1 } space by vector fields that are divergence free, i.e., locally 
we need that (f) 1 — > (f) 1 + £ J with Dj £ 7 = 0. Here T>j denotes the covariant derivative 
in the manifold parameterized by {(p 1 } which we call M.§. The symmetry we demand 
is invariance under volume preserving diffeomorphisms of M.<j,; this symmetry group 
is often denoted as Sdiff(TW^). 

Before proceeding to construct an action with the symmetries described above 
let us note one consequence of the field reparametetrization invariance. Consider the 
current one- from obtained by taking the spacetime dual of the volume form of M.^: 

J = * (dct) 1 A dcj) 2 A • • • d^ 1 ) , (2.2) 

or directly in components: 

d-l 

J" = y—fy ^- ad - 1 ^ 1 -/ d - 1 J! ^ ■ (2-3) 

This current is is trivially conserved 

V^ = 0. (2.4) 

We will soon see that this current is to be viewed as the entropy current Jg and the 
volume preserving diffeomorphism symmetry ensures this conservation automatically. 

We will find it convenient to split the current J M into a scalar and a normalized 
(i-vector: 

J^ = su^, s = yf-J^ , u^u^ = -l (2.5) 

and build the Lagrangian out of these fields (respecting spacetime diffeomorphism 
invariance). The notation is intentionally suggestive: s is to be identified with the 
entropy density of the fluid and with its normalized d-velocity. 

2.2 Operator dimensions: setting up the gradient expansion 

Now we are left with the task of writing down the general effective action with the 
symmetries of the fields. To get off the ground, we need to understand the canonical 
operator dimensions for such a construction. Since the (p 1 fields are to be viewed as 
Goldstone modes for the embedding of the fluid into the background spacetime, they 
turn out to have mass dimension = — 1. This effectively implies that we treat 
them as phase fields and consequently [dcp 1 ^ = implying that J M is a dimensionless 
operator in the field space (and a tensor density in physical spacetime). 

Our task is now clear: we should first classify all operators Oa with canonical 
scaling dimension A = 0, 1,2, . . . etc.. However, since [<i0 7 ] = even an operator 
with A = can have arbitrary many derivatives when we attempt to build it out 
of the fundamental fields 7 . Fortunately, the invariance under SdiS(M^) comes to 
our rescue. We use this large symmetry to make the following assertion: 
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Claim: The action S is a functional of the entropy current and its derivatives, 
the latter being treated as the zeroth order operator. The full action can be written 
in derivative expansion as an infinite sum 



where i th order the action Si contains % space-time derivatives. In addition to the 
operator J 11 (equivalently s and u^) we should also keep track of terms that involve 
operators built out of the background geometry. In the present discussion of neutral 
fluids these are simply operators built out of the intrinsic geometry of the background, 
and hence involve the metric tensor g^ v and the curvature tensors built from it. 

We now proceed to analyze the effective action in detail at the first three orders in 
the spacetime derivative expansion, based on the above assertion without providing 
a definitive proof of it. Note that one really should classify all invariants under 
Sdiff(Al^) at a given order in the derivative counting described above; it is the latter 
that ensures that the dynamics of (f) 1 is equivalent to energy-momentum conservation 
while guaranteeing entropy conservation. 

2.3 Zeroth order: Ideal fluid hydrodynamics 

At zero derivative order the only scalar that can be constructed out cf) 1 respecting 
Sdifl^A^) invariance is the norm of J M which according to (2.5) is proportional to 
the square of the entropy density. Therefore the zeroth order action can be written 
as (the normalization is to ensure a simple stress tensor) 

So = -2 / J=g f(a) (2.6) 



where f(s) is some arbitrary function of the entropy density. 

Given this effective action, we can vary with respect to the background spacetime 
metric using (A. 2) to obtain the stress tensor 3 



= (sf'(s)-f(s))g^ + sf(s)u»u» (2 ' 7) 

In the last line we have identified the energy density e(s) = f(s) and used thermody- 
namics to identify p(s) = (s f'(s) — f(s)) with the pressure. P^ IV = g^ v + u fl u u is the 
projector 4 in the direction transverse to u M . After this identification the stress tensor 



Henceforth for any function F(s) we denote by F' the derivative with respect to entropy density 
s Le, F>(s) = f . 

4 Our conventions for fluid dynamical tensors are collected in (3.2). 
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presented in (2.7) has taken the form of the ideal stress tensor for an uncharged fluid. 
Thus the Lagrangian formulation for an ideal fluid simply involves writing down the 
energy in terms of the entropy density. 

Let us quickly check the equations of motion: extremizing this action with respect 
to the (ft 1 we get: 



#a / 
{u v P^-u a P^)V a [f{s)u v } = 



{21 



In the second line we have used the fact that J M V M 0/ = w M V M 0j = 0. 

There are thus d — 1 Euler-Lagrange equations (one for each (p 1 ) which we can 
conveniently assemble into the conservation equations of the energy-momentum ten- 
sor projected transverse to the velocity The equations in (2.8) automatically 
follow from simply P MCf V U T^ = 0. Further, as already argued in [16] the remain- 
ing conservation equation u^V^T^ = is implied by the conservation of entropy 
current J^ 1 and therefore identically zero. Therefore the equations of motion for the 
basic <pi fields have no more content that that of the conservation of stress tensor as 
desired. This is of course in keeping with the general motivation for working with 
the fields (j) 1 ; the only novelty is that Sdiff(jM^) invariance allows us to trade energy 
conservation for entropy conservation. 

2.4 First order: absence of viscous effects 

At first order in derivative expansion we can construct one non-trivial scalar out of 
J 74 , namely J^V^s. 5 Therefore naively we can write an action at first order 

Si = J V=gfi(s)J^,s = J V=gJ»VMs) (2.9) 

with for some = fi(s). This term is however a total derivative, for we can 

rewrite Si as 



Si 



f [V„ r) - h{s) V»J»] (2.10) 



with the second term being identically zero. 

In our analysis we shall ignore all the total derivative terms in the action, since 
they don't contribute to equations of motion. It then follows that Si can be ignored 
in what follows; therefore the first order correction to the stress tensor is also zero. 
The first non-trivial correction appears at second order. 6 

5 The one other scalar built out of J is its divergence of J p which is identically vanishes. 

6 This statement is not true in three dimensional fluids with parity violation: we can construct 
a one-derivative operator hi(s) e^ vp J^ V v J p which is not a total derivative. This contributes to 
transport; see [15] and also Appendix B for details. 
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This is of course intuitive; we are describing a non-dissipative fluid with an 
exactly conserved entropy current. At first order we know that we can have vis- 
cous contributions to the stress tensor, but these generically produce entropy, since 
V M Jg oc 77 a^ u + ( 2 where a^ u and 6 are the fluid shear and expansion respec- 
tively, cf., (3.2). Requiring entropy conservation for arbitrary configurations forces 
upon us 77 = ( = 0. 

This point is likely to make the reader uncomfortable; while we know that fluids 
with ( = are physical (e.g., conformal fluids), vanishing of 77 seems a bit strange. 
We postpone the discussion of whether this is physically acceptable to §5. 

2.5 Second order: novel non-dissipative fluids 

We now turn to the second order in gradients and argue that we have a five distinct 
on-shell inequivalent operators that we can write down in the effective action. First, 
note that now we encounter two types of scalar operators: ones that involve explicit 
curvature tensor and the others that do not. 
There are two curvature dependent scalars: 

R and R pv u ti u u , (2.11) 

which can enter our effective action multiplied by arbitrary functions of s. 

The scalars without curvature can again be divided into two subclasses, ones 
where both the derivatives are acting on the same J M (e.g. J M V 2 J^) and the others 
where the two derivatives are acting on two different J M (e.g. V M J V V^ J u ). But any 
term which is of the first type can always be recast into a term of second type upto a 
total derivative. Therefore for the second order action we only need to classify those 
scalars which are constructed as product of terms with single derivative acting 

Equivalently we can view the construction in terms of s and and ask for 
scalars that are built out of the one-derivative vector V M s and the one- derivative two 
tensor V^tv A-priori we would write down six such scalars: 

(V^XW) (2.12) 

where for any vector {A^) 2 denotes (A^A^) for brevity. We now argue that in 
fact from the basis of 8 operators contained in (2.11) and (2.12) we only need to 
consider 5 operators. 

Generically the action constructed at second order will have a field redefinition 
ambiguity. In order to fix this we can work only with those operators which are 
on-shell inequivalent. The equations of motion arising from the zeroth order action 
So (2.6), together with V^J^ = allows us to relate derivatives acting on s to 
derivatives acting on u^. So we only need to consider two-derivative scalars obtained 
from contractions the two-tensor V^m^. Using the decomposition of this two tensor 



(n a V aM „)V^, (V„0(W), 
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into a transverse vector (acceleration), symmetric traceless & anti-symmetric two- 
tensors (shear & vorticity) and a scalar (expansion) (cf., (3.2)) i.e., 



d-1 ""• 



we find a basis of four scalar operators: 



„2 Jl , ,2 

a , a , w 



2 . 



(2.13) 



(2.14) 



Furthermore, using the the fact that R^ v u^u u is given by a commutator of derivatives 
acting on w M u v , we can eliminate it in favor of operators built from the velocity 
(2.14), see Appendix A for details. This is physically the most intuitive basis given 
the observations earlier about rj = ( = 0. Hence the action at second order finally 
has five terms. 

From a computational perspective however, we find it convenient to switch to a 
different basis which is more efficient for variation of the action with respect to the 
metric. Introduce thus the second order action: 



S 2 



-g 



K^s) P afi {V a u u ) {V p u u ) + K 2 {s) {V„u v ) (Vtt") 
+ K 3 (s) (J^V^s) 2 + K 4 (s) V M s V^s + K 5 (s) R 



(2.15) 



which satisfies all our physical requirements. 

Varying this action with respect to the metric using (A. 2) we get stress tensor 
contribution: 



IF" 



1 SS * _ yr^U 

r=9 Sg„ ~ ^ U W 



(2.16) 



8=1 



where 



u (l) 



K x - sK[ 



P aP {V a u e ){Vpu 6 



pfiv _|_ 



+ tfiaVV - K 1 (V M M a )(V I/ u Q ) + KiP^V^XV^) 

(2.17a) 



U (2) 



K 2 - sK' 



(V ««)(V^«' 



(V^XW) + ^V" (^ 2 V Q ^) 



(2.17b) 



-TTflV 

U (3) 



3^3 + , _ 



pM" _|_ 



S 4 A ? 



e 2 



u"u" (2.17c) 
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U (4) 



JTflV 

U (5) 



K & + sK' 



(Vs) 2 + sK 4 V 2 s 



plLV _ 


[(f) H 







R 



K, - sKL 



2 



(2.17d) 



(2.17e) 

We note that some of the expressions above can be simplified upon using the standard 
decomposition of the covariant derivative of the velocity field u^; see Appendix A. We 
have also used the entropy conservation equation in fl^ to express J M V M s = — s 2 0. 



That total stress tensor T^ v = + U^ u is of course conserved and by our 
earlier arguments this conservation implies the Euler-Lagrange equations for the 
action So + 5*2. Thus using this effective action approach we have constructed a five 
parameter (in the sense described in footnote 1) family of non-dissipative fluids. 

3 Constraints arising from a conserved entropy current 

We now turn to a different approach to describing non-dissipative fluid dynamics 
based on an analysis of the entropy current. This approach has been used to investi- 
gate neutral fluids initially in [3] and more recently [7] has analyzed the constraints 
on regular dissipative fluids exhaustively up to second order. In this section we shall 
see how the existence of an exactly conserved entropy current constrains the possible 
form of the stress tensor at first and second order in a gradient expansion; here we 
shall follow analysis of [7]. 7 

To begin with we have to use symmetry and on-shell equivalence to have a naive 
count of the number of possible terms in a hydrodynamical stress tensor at first and 
second order in gradients. This has been done in [7]; it turns out that the stress 
tensor at first order has 2 transport coefficients and at second order has 15 transport 
coefficients. Writing T^ v = + W J,V the most general form of IF 1 ' up to second 
order in gradients (constrained only by symmetry and on-shell equivalence) is given 
by the following expression: 



n 



+ T 



T U a V a (T fa-,,) + K\ Rfav) + 1^2 F^ v ) + A © CT M z, 

+ Ai (T( M Q a av ) + A 2 cr ^ a to au ) + A 3 to av ) + A 4 a^a^ 



(3.1; 



+ T P. 



Cm a V a Q + ( 2 R + ( 3 R 00 + £i© 2 + + £ 3 ^ 2 + £ 4 a 2 



7 While the analysis described here can be generalised to arbitrary dimensions, we will restrict 
attention to d = 4 in this section mainly to avoid accidental relations (especially in d = 2, 3) which 
affect our choice independent tensor structures. 
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with T being the (local) temperature. In writing this we have introduced some 
notation which are defined as follows in d — 4 space-time dimensions: 

= The normalised four velocity of the fluid 
pv-v = _|_ u n u v _ p ro j ec tor perpendicular to 
= V a u a = Expansion, a u = u a V a u a = Acceleration 



a , V = p.apu, ( ^aU, + V,U a _ jL_g a ^ = Shear tensor 



(3.2) 



a 2 = a au a^, to 2 = w^w"" 



with being the Riemann tensor of the background geometry. Furthermore, we 

define a projection of any tensor A uv onto its symmetric transverse (to w M ) traceless 
part via 

A a /3 + A/3 a 



A ab P ab 
d- 1 



o , , 9a() ) • (3-3) 

A conserved entropy current is the special case of the entropy current with non- 
negative divergence. The constraints on the second order transport coefficients due 
to the existence of a positive divergence entropy current have been determined in 
[7]. The calculation in this section involves a slight modification of such a general 
analysis (which is of course valid for all fluids). In fact the analysis up to Section 4 
of [7] can be adapted unchanged. Subsequent discussions will differ, but in a very 
simple way, so that the results can easily be read off from the equations presented in 
Section 5.1 and Appendices A and B of [7]. 

Let us first briefly recall the logic used in [7] to determine the constraints on 
transport. We shall then be able to see how to modify it to avoid dissipation and 
obtain the constraints on the transport coefficients. The upshot of our discussion 
will be that we will show that of the {2 + 15} transport coefficients appearing in 
(3.1) at first and second order respectively, one fixes {2 + 13} in terms of the five 
arbitrary functions appearing in the entropy current. 8 

Entropy current has two parts, one is a canonical piece and the other is the 
correction. The canonical piece is completely fixed in terms of the zeroth order 
entropy current and IP y , the gradient correction to the ideal stress tensor. 

u v IP 1 ' 

= Jean + J cor J ^can = S U ^ p (3-4) 

In [7] a particular fluid frame (Landau frame) was been chosen by demanding 

u u U^ = (3.5) 



8 We will denote the number of arbitrary functions appearing in various quantities of interest at 
first and second order in gradients by ordor + # 2 ,,ci order} f° r simplicity. 
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In this frame J^ an is equal to entropy density times the local velocity 

The correction to the zeroth order canonical entropy current is determined using 
the fact that its divergence should never go negative on those fluid flows which satisfy 
the equations of motion. This argument remains unaffected even if we demand the 
total divergence of J M to be zero. 9 Therefore we have the same Jjf or as in Eq.(1.4) 
of [7]. The final form of <P in Landau gauge (correct up to second order) which we 
should consider is given by: 

J" = s u» + V„ [2 Ax u [fl V" ] T] + V U {A 2 T uT) 

+ A 3 (bT - \g^R\ u u + + [e WT - P a ?VpU» V a T] (3.6) 

+ (Bi u 2 + B 2 Q 2 + B 3 a 2 ) + B 4 [V a s V a s + 2 s V^s] 

Given such an entropy current one then computes its divergence. The latter also 
is best viewed in the decomposition of the canonic contribution and the correction 
piece. Using the equations of motion for + IP" together with thermodynamic 
relations one finds (in d — 4): 

= v M J£ n + v^ or = -i (a^ ur + 1 p m + v,j» or (3.7) 

The key point to note is that the divergence of the canonical piece always involves IP" 
as one of the factor and therefore will contain the information about the transport 
coefficients. 

For generic fluids in [7] the idea was to rewrite (3.7) as a sum of perfect squares 
(up to fourth order in gradients) so as to ensure positivity. All the terms which 
could not be cast into the perfect square form were set to zero - this gave the final 
constraints on the transport coefficients. 

However, here we are interested in a zero divergence entropy current; this is the 
place where the present computation differs from that of [7]. We now proceed to 
outline the consequences of demanding that V^J M = 0. 

Note that the two terms in (3.6) with coefficients A\ and A 2 are explicitly of zero 
divergence and drop out of our analysis. The rest of the five terms in J^ or with five 
independent coefficients have non-zero divergence and therefore have to be cancelled 
against corresponding pieces coming from the divergence of J^ an and constrain IP". 

When we compute the divergence of the entropy current given in (3.6) up to third 
order in gradients and express the final answer in terms of the on-shell inequivalent 
data, each term turns out to have either a^ v or as a factor. One can thus argue 

9 A-priori the most general entropy current for parity invariant fluids up to second order in 
gradients has {2 + 13} parameters. By examining the number of independent scalars produced 
when one considers V M J p we find that {2 + 6} of these have to be set to zero leaving behind the 
{0 + 7} parameter set quoted in (3.6). 
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that schematically (please consult Appendix B of [7] for explicit expressions): 



V A4 J£ r = <v# u ' + efl (3.8) 

The final expression of the divergence of the full entropy current is given by the 
following. 

V = -<v L 4--B^)-Q - B (3.9) 



x T J V 3T 

Generically the RHS of (3.9) will be zero if we set 

n (M = rp B (^) + q^u and n ^P^ = 3 T B (3. 10) 

where Q^ v is any symmetric traceless tensor satisfying 

<y^Q' w = o 

Using the schematic equations (3.7), (3.9) and (3.10) we can draw a few imme- 
diate conclusions. 

• From (3.10) upon demanding zero- divergence to third order, naively it appears 
that we can determine the traceless part and the trace part of the stress tensor 
separately in terms of the parameters appearing in the entropy current. This is 
generically true but with an important exception: when both the traceless and 
the trace part of the stress tensor contribute the same term to the divergence 
we can only fix a linear combination. This occurs for the terms Ao 6 0" M „ in 
the traceless part and £ 2 cr 2 in the trace part of our stress tensor (3.1). Both 
of these end up contributing a term proportional to a 2 in the divergence of 
the canonical entropy current which is to be cancelled by the divergence of the 
'correction' part; this is manifest by writing: 

A e <v + 6 P,u o 2 = (6 <v + P, u a 2 ) + (e ^ - P, u a 2 ) 

(3.11) 

Therefore our analysis will only be able to constrain Ao + £2- 

• Since J^ or does not contain any term which is of first order in derivatives, both 
the tensor B^ v and the scalar B contain exactly two space-time derivatives. 
It then immediately follows from (3.10) that the first order correction to the 
stress tensor has to be zero, i.e., both shear viscosity 77 and bulk viscosity £ 
have to be zero, as is physically sensible for a non-dissipative fluid. 

• For parity invariant fluids Q^ u is proportional to a^ a uj av ^ at second order. 
Hence from (3.10) it also follows that the transport coefficient A2, multiplying 
this term in IF^, will remain unconstrained by the entropy current analysis. 
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We therefore can immediately conclude that {2 + 13} transport coefficients ap- 
pearing in (3.1) are fixed in terms of {0 + 5} non-trivial arbitrary parameters ap- 
pearing in J[? or . The explicit relations are given by: 



dT 

A + 6 = ( B 3 - s 
dB 5 



2B 3 , «! = A 3 , 
dB 3 



k 2 = T 



djh 
dT 



ds 



ds 



Ai 
A 4 
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dT 



A, = T 



dEk 
dT 



4B 1 
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d 2 B 5 
dT 2 



dB. 



+ T ^ + 2 B A T 2 . 

dT \dT 



ds 



(3.12a) 
(3.12b) 
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(3.12d) 
(3.12e) 
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(3.12g) 



where ^ = T , dT 

entropy current. To summarise have found a {0 + 7} parameter family of such fluids; 
5 parameters appear naturally in our parameterisation of the entropy current and 
two of the parameters that appear explicitly in the stress tensor viz., A2 and a linear 
combination of Ao,^2- We note in passing that Ai = k 2 (which in turn in related to 
Ki) for absence of entropy production. 

So far our discussion has been restricted to generic neutral fluids. One can fur- 
ther constrain the class of hydrodynamic theories by demanding that the fluid be 
conformal. Since the trace of the energy momentum tensor has to vanish for confor- 
mal fluids, this in particular implies that the only non-trivial transport coefficients 
are (cf., [7, 19]): 



T 



3 An 



K,2 = 2 K\ 



Ai 



A 5 



A, 



(3.13) 
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All other second-order transport coefficients must vanish. A-priori since we have the 
remaining ten transport parameters determined in terms of five arbitrary functions 
one might suspect that there is no non-trivial solution. Working things out however 
we find that choosing 10 

A 3 = a 3 T } B 1 = b 1 T, B 3 = b 3 T A 2 = £ 2 T, (3.14) 

ensures that all the conditions from = are satisfied. The rest of the parameters 
B 2 and _B 4 are determined in terms of a 3 , in particular, £> 4 = — ^5 and B 2 = — | a 3 T. 
Note that scale invariance dictates that s oc T 3 . Changing basis of variables, the non- 
trivial transport coefficients can be taken to be r, Ai, A 2 and A3 (since k is given 
by Ai by (3.12)). We thus claim that the entropy analysis leads us to conclude the 
existence of a four parameter family of conformal fluids. 11 

4 Frame invariant formulation and comparison 

We now turn to a comparison of the analysis of the preceding two sections of non- 
dissipative fluids. Whereas in §3 we examined the constraints arising form the ex- 
istence of the most general zero divergence entropy current, the analysis of §2 was 
predicated upon the existence of such a conserved entropy current. Therefore the 
stress tensor derived in §2 must be a special case of the stress tensor derived in §3 
up to a identification of parameters (we have conveniently recast the effective action 
formalism in terms of the fluid variables already); the same must also hold for the 
entropy current. 

However, the analysis in the previous section has been done in Landau frame 
(3.5) which defines the velocity w M (and temperature). Unfortunately the stress 
tensor obtained in §2 by varying the second order correction to the action is not in a 
Landau frame (this can be explicitly checked using the explicit form of IP 1 ' given in 
(2.16)). In fact, we can view the analysis of §2 to be in the 'entropy frame' defined 
by the condition that entropy current does not receive any gradient correction, i.e., 

= s/ , at all orders (4.1) 

Given these observations one should not directly compare the given in (2.16) with 
the expression derived in §3. We have to perform a field redefinition, i.e., redefine 
the velocity and temperature of one of the expressions so as to bring both answers 
into a common frame. 

10 These relations can be read off directly from the analysis in Section 6 of [7], cf., Eq (6.7) of 
that paper. Noting that scale invariance also fixes Ai = 2 03 T and A2 — 0,2 T, the only difference 
from the analysis of [7] is the constraint on B^: we find here that B2 = — | 03 T, which differs by a 
factor of —2 from the value found earlier (for comparison please note that a' horc — 2a 3 lere ). 

11 Since scale invariance fixes the temperature dependence of our transport coefficients, we really 
have only four parameters for conformal fluids. 
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Happily, one can avoid doing an explicit field redefinition by choosing to com- 
pare only those combinations of IT^ which are invariant under such redefinitions of 
velocity and temperature (i.e., a frame transformation). Moreover the frame invari- 
ant combinations for both charged and uncharged fluids and superfluids have been 
worked out in [5]. Their formulation is applicable at the first non-trivial order of 
corrections to the stress tensors (and currents) in gradient expansion. Generically 
it is applied to the first order piece of IP" as we encounter viscous effects at that 
order. But we have already seen that for uncharged non-dissipative fluid the first 
non trivial correction appears at second order, which is the one we are interested 
in. Therefore in this particular case we can use the frame-invariant formulation to 
simplify the computation despite working at second order. 

For uncharged fluids the frame invariant combinations that contain the infor- 
mation about the transport coefficients are a transverse traceless two tensor and a 
scalar: 

Of = P^P^U^ - - P^ P a/3 Ii aP , 

_ P^IV s (dT\ ( 4 ' 2 ) 



3 T \ds J 



We compute Cf and C 2 using the stress tensors computed in §2 and §3 and compare 
the two answers. 

Given that stress tensor derived in section 3 contains 7 arbitrary functions and 
that derived in (2.16) has 5 functions, we expect the frame invariant combinations 
to match provided we express the 7 transport parameters {A 3 , B±, B 2 , B 3 , B±, \q, \ 2 } 
in terms of the parameters ifj (for % = 1, . . . , 5) appearing in our effective action 
analysis of §2. A somewhat tedious algebra 12 reveals that 

A B = K2 ~ Kl B = Kl + K2 

K x + K 2 s 2 dT dK 5 s 4 T ^ 

Do — + — A 1 ! 

2 6T T 2 ds ds 2T 3 



K 4 1 dT dK 5 



2T T 2 ds ds 



T ds 3 dT T \ ds ds 

A^2^1 (4.3) 

Curiously A2 which does not enter into the entropy current analysis is very simply 
related to B\\ from the above A2 = — 4£> 3 for non-dissipative fluids arising from an 
effective action. 



12 We need to use the zeroth order equations of motion in (2.17) to extract the correct tensor 
structures; in particular we make use use of V M s = sQu^ — T a M . 
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In a generic frame the expression for the full entropy current is given in (3.4). In 
fact in this expression one can show that — ^u u IF 1 ' is frame invariant, thereby 
demanding that J^ or is insensitive to field redefinitions. Given that the analysis of 
§2 is in the entropy frame, one is led to the following conclusion: ^ u v IP" computed 
using (2.16) should be equal to J^ or of (3.6) once we substitute the identifications as 
given in (4.3). This algebraic check works out as expected and along the way we get 
a useful bonus. One is now able to fix the two remaining parameters of the entropy 
current A\ and A 2 which by virtue their vanishing divergence remain unconstrained. 
One finds 

1 dK b K 2 -K l . . 

Ai = - T -^. * = -V^ (4 ' 4) 

Notice again the curious fact that A 2 which also involves the vorticity is related to A2 
in a very simple manner T A 2 = 2 A2 = \ B\. It would be interesting to understand 
why these relations are forced upon the fluids that arise out of the effective action 
constructed using Sdiff(A^0) invariance. 

Finally, note that we can fix the functions Ki to ensure that the fluid is conformal, 
by using the constraints on the functions entering the entropy current analysis and 
(4.3). We find: 13 

K x = -(6i + 6 3 )s5, K 2 = (61 + 63) si , 

2 1 2 2 

K 3 = --b 3 s-^ , K 4 = --a 3 s-^, K 5 = -a 3 s 2 z, (4.5) 

accounts for all the constraints. Since now l 2 which determines A2 is fixed in terms 
of 63 from (4.3), we learn that the conformal fluids which arise from an action are 
required to satisfy one additional relation. Using the parameterization (3.14) the 
transport coefficients for conformal fluids arising from the action principle can be 
encapsulated as: 

t = 3 A = (a 3 + 63) T , k 2 = 2 ki = Ai = 2 a 3 T , 

A 2 = -46 3 T, A 3 = 2(a 3 -2b 1 )T (4.6) 

Thus, the action formalism only allows us to explore a three parameter family of 
non-dissipative conformal fluids. 



13 Notice that we are not a-priori demanding scale invariance of our action, but rather using the 
identification between the action and entropy current analyses (4.3) to constrain the functions Ki 
appearing in the action (2.15). It is easy to convince oneself that the counting works out correctly 
by noting that the conformally covariant scalars are three in number: these are built out of cr M „ <7 M!y 
uj^ u uj ij,v and a specific linear combination of {8 2 ,o 2 ,i?}. We have refrained from translating this 
into the basis we chose to work with preventing a direct comparison immediately; the constraint on 
Ki (4.5), should follow directly from this translation. 
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5 Conclusion 



Motivated by a need to understand the fundamental constraints on an autonomous 
theory of hydrodynamics, we have in this note explore neutral non-dissipative fluids. 
In particular, we have examined the conditions on such fluids (neutral) arising from 
demanding the existence of a divergence-free entropy current as well as from an 
action principle. While both these analyses are consistent with each other, we find 
a larger parameter family of non-dissipative fluids from an entropy current analysis. 

The entropy current formalism determines {2 + 13} out of a total {2 + 15} 
transport coefficients of the first and second order stress tensor in terms of the 5 free 
parameters appearing in the most general zero divergence entropy current (thereby 
predicting 8 linear relations among 13 of second order transport coefficients). Two 
second order transport coefficients do not enter into the entropy current analysis: 
these are A2 which multiplies the (shear) x (vorticity) contribution and a linear 
combination Ao — £2 which multiplies a trace-free combination of expansion and 
shear. All in all we have a seven parameter family of second order non-dissipative 
fluids which are consistent with absence of entropy production for arbitrary flows. 

On the other hand, in the effective action action approach the existence of a 
zero divergence entropy current is ensured by demanding the reparameterization 
invariance in field space Sdiff^A^). We have found a five parameter family of effective 
actions at second order, which lead thence to a stress tensor where all the {2 + 15} 
transport coefficients are determined in terms of these 5 parameters. The resulting 
stress tensor is of course a special case of the one predicted by the entropy analysis 
and upon suitable identification of variables one is also led to fixing the free transport 
coefficients A 2 and A — £2 of the latter construction. 

The main outstanding question of our analysis is why the effective action ap- 
proach gives a smaller parameter family of second order non-dissipative fluids than 
that predicted by demanding the presence of a zero divergence entropy current. A- 
priori one can think of two distinct reasons for this mismatch: 

(i) We have not identified all possible dimension two scalar operators built out of 
the / , the basic variables entering into the effective action. As we discussed 
in §2 we worked out the independent structures that preserved the volume 
form on the (p 1 manifold Ai^ using structures built out of the current J M and 
its derivatives. It is plausible (though we think unlikely) that there are other 
admissible structures. 

(ii) More intriguing (assuming we have identified all possible contributions to the 
action) is the possibility that there are indeed further constraints on fluid dy- 
namics coming from the existence of an effective action, which are not simply 
captured in terms of symmetry constraints on the stress tensor and demand- 
ing the local form of the second law (or rather zero divergence for the case of 
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non-dissipative fluids). If this were the case, then we would have an interesting 
window of opportunity of using these non-dissipative fluids to learn about the 
potential constraints on hydrodynamic expansions. Perhaps out results hint at 
additional microscopic constraints on the transport coefficients A2 and Ao — £2? 

The mismatch in the parameter count is the salient feature of our analysis and 
its shows up at second order for neutral fluids. We would be remiss to not point 
out that the effective action approach can be seen to give a restricted class of fluids 
when we allow for parity-violation already at first order. The authors of [15] have 
previously used this approach to describe Hall viscosity (a new transport coefficient) 
in three dimensions. This parameter multiplies the parity-odd e ap ^ p u a Op term in the 
stress tensor. One can show that the coefficient here is entropy preserving. In fact it 
does not enter into the entropy current analysis (for pretty much the same reason as 
(Ta* bj v ^ a in our analysis of §3). Curiously this term is again not reproduced from the 
canonical local action one writes down. Moreover, the entropy current following the 
recent analysis of [20] (who also included a conserved global U(l) current) reveals a 
two parameter family of parity-odd dissipation less fluids, while the effective action 
predicts a single parameter (see Appendix B for details). Likewise in the analysis 
of [17] who explored the possibility of describing anomalous charged fluid transport 
in d = 2 using this effective action approach, one obtains a one-parameter family 
of anomalous charged fluids at first order in gradients. On the other hand entropy 
considerations [6, 21] as well as more recent analysis using equilibrium partition 
functions [22] clearly show a two parameter family of such fluids (only one of the 
parameters is related to the global current anomaly). 

It would be interesting to extend this analysis to charged fluids, superfluids and 
fluids with anomaly (in higher dimensions) to see whether we obtain more constraints 
from demanding an effective action's existence. 

At the end of the day if we believe the implications of the effective action for- 
malism, we should be able to derive the constraints on transport directly from mi- 
croscopic theory. A useful avenue for exploration is perhaps analysis of stress tensor 
correlation functions. For instance as shown in [23] (see also [24] for non-conformal 
fluids) one can derive Kubo formulae for the second order transport coefficients and 
in particular for A2. One should examine whether the condition of vanishing 7] forces 
constraints on the particular three-point function determining A2. 

While still on the topic of fluid dynamics, in [25, 26] the authors have written a 
partition function or a generating functional for stress tensor evaluated on the most 
general time independent fluid flow on any arbitrary but static background. One can 
check using the action (2.15), that there is a three-paramater family of equilibrium 
actions (essentially a^ v = = for stationarity), which agrees with the counting 
based on sources that be turned on maintaining time independence. It would be 
interesting to explore the connection between these two formalisms. In particular it 
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would be nice to know how one can determine the partition function describing the 
time independent situation from the action written for a time dependent but strictly 
non-dissipative fluid flows. 

The bulk of this paper has been concerned with analysis of fluid dynamics in 
its own right as an effective action. We have not so far touched upon the issue of 
whether non-dissipative fluids are physical, nor have we mentioned any connections 
to holography. Now is the time to remedy these lacunae in our discussion. 

As we remarked in §1 the fact that on-dissipative fluids demand t] = ( = makes 
one suspect that this type of hydrodynamics cannot arise from a sensible unitary 
quantum field theory. Indeed while ( = is realized in conformal fluids, one expects 
on general grounds the viscosities to be bounded from below from general arguments. 
For instance using the uncertainty principle one can generically argue that i] > a s 
where for a change we have explicitly indicated the fundamental constants. There 
is an undetermined constant a ~ 0(1) in the above analysis. Indeed it has been 
speculated following the seminal work of [27] that there is a lower bound on shear 
viscosity; the status of this bound has been the subject of much scrutiny and the 
current wisdom is that 2 > where 7 is a fraction which seems to depend on the 
details of microscopic theory (see [28] for a recent review of the bound). Likewise, 
while ( = is forced upon one by scale invariance, in non-scale invariance theories 
it is also conjectured that there is a lower bound on (, viz., ( > 1r\ (^-j- — c^) with 
c s being the speed of sound [29] (this bound is much less explored). Based on these 
arguments one might indeed be tempted to take the viewpoint that non-dissipative 
fluids of the type discussed herein are unphysical. This is supported for instance 
by holographic computations of shear viscosity in higher curvature gravity theories; 
[30, 31] studied rj/s in Gauss-Bonet-AdS gravity and argued that while one could use 
the Gauss-Bonet coupling to lower 77 (in fact all the way to zero) demanding sensible 
causal properties of the bulk (most likely related to unitarity in the field theory) 
results in a lower bound for rj away from zero. It would be interesting to explore this 
example in some detail to understand whether the fine tuned Gauss-Bonet theory 
with rj = has non-trivial second order transport coefficients which refrain from 
entropy production. 

Finally, one of the reasons to take this explorations seriously is the potential it 
has to teach us about honest dissipative fluids. One might ask how could one deform 
away from the situation where entropy is conserved to one where the local form of 
second law is valid in standard form. A clue comes from the origin of the identity 
J M = 0; this is forced upon us by the fact that we demand invariance under volume 
preserving diffeomorphisms on <ft space. Suppose we were to relax this condition to 
just demanding diffeomorphism invariance - this is a sensible requirement to impose 
since we are still allowed to relabel fluid elements arbitrarily. Can one then write 
down an effective action for dissipative fluids? The tentative answer coming from the 
analysis of [32] is yes: by using holographic intuition for the renormalization group 
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these authors have constructed a leading order effective action which seems to capture 
the effects of shear viscosity. Indeed as they point out in the holographic context 
the fields (p 1 which we view abstractly as labels for the fluid, pick up a fascinating 
geometric meaning. Since the fluid/gravity correspondence [33] asserts that arbitrary 
fluid configurations of a field theory with a holographic dual is given by a black hole 
solution in an asymptotically AdS spacetime with a regular event horizon, 14 we can 
parameterize our solution in terms of the spatial geometry of the horizon which is 
then Lie transported along the future horizon generator. The natural coordinates on 
these spatial sections are nothing but the 7 : in equilibrium (or indeed in the absence 
of dissipation) this follows from the analysis of gravitational entropy current of [35] . 
It would be interesting to explore these relations in greater detail and we hope that 
future studies enables one chart out a clear autonomous theory of hydrodynamics. 
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A Some useful results 

Conventions: We collect some useful formulae in this appendix. First with regard 
to conventions: in the text we use the standard (anti-)symmetrisation convention 
i.e., = | (A pu + A up ). We define the curvature tensors as: 

R p a/3u = dpY p ai/ — d u Y p a p + r^r^g — r^r^, , r^ v = r p ^ (a.i) 

14 This assertion extends beyond AdS to encompass long-wavelength world-volume fluctuations 
of black branes as exemplified in the blackfold approach [34] . 
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Useful variational formulae: The reader will find the following useful to derive 
(2.16): the variation of J M , u^, s and the connection) with respect to the 

metric are given by 



5J» 

5s 
STL 



f ) 4* <W 



u a vP 5g a/3 



s 
2 



(A.2) 



P afi 6g a p 



Eliminating the velocity projected curvature scalar: To show that R^ v u u 
can be eliminated on-shell note that 

J Jg K(s)R^u v = -jV9 K(s) [V„, V p ] u»u" 
^ K{s) [u v V p V v u p + u a V a (u^Vps)] 

(A.3) 

Vg{if( S )(V/)(Vy) + [(u.V)u"]d p K(s) + (u.ds)V^ [K(s)u»]} 
V9^M(V P /)(Vy)-T^[( M .VKf + ^ (u.ds) 2 } 

Simplifications of the stress tensor derived in §2: Some of the expressions 
appearing in the stress tensor derived from the action can be simplified when we use 
the decomposition (2.13) and also the zeroth order equations of motion. The latter 
can be cast in a useful form which reads: 



ds 



(A.4) 



For reference we record that the combinations appearing in (2.17a) and (2.17b) 
can be simplified into the form: 
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(A.6) 



This is sufficient to read off some of the transport coefficients, especially r directly. 
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B Parity-odd neutral non-dissipative fluids in d = 3 



In this appendix we will consider parity-odd neutral fluids in 2 + 1 dimensions at first 
order in derivative expansion. As we have seen in §2 there are no non-dissipative 
terms at first order in hydrodynamics in 3 + 1 dimensions. This fact is also true in 
2+1 dimension, if we preserve parity. However if we consider parity violating effects 
then it is possible to write a first order term in the action using our basic variables 
(j> , which is manifestly invariant under reparameterizations of the the fields [15]. The 
desired term has the form 



^(parity-odd) = J d 3 x ^^ ^A^ VvUx = j d 3 x ^^ ^ (B.l) 

where <r(s) is some arbitrary function of the entropy density which will ultimately 
determine the transport coefficients appearing in the stress tensor. Here 

Q = t puX u^ v u x (B.2) 

is the vorticity which is a scalar in 2 + 1 dimensions. Our conventions are that 
e puX = ^= e puX where e pup is the flat space Levi-Civita tensor, with orientation 

The first order corrections to the stress tensor that follows from the action (B.l) 
is given by 

(parity-odd) 2 (B 3) 

+ 2 q(s) e^ px u^V p u x + q'(s) e^ px u x V p s u v \ 

Note that there is no frame invariant genuine tensor in this expression of the stress 
tensor. In other words if we take ^p arity _ odd ) from (B.3) and evaluate using (4.2) 
suitably modified for 2 + 1 dimensions, we find that it evaluates to zero. There is 
however, a frame invariant scalar data in (B.3), since for this correction to stress 
tensor C 2 in (4.2) evaluates to 

ct arity " odd) = - (^00 + ~ ?(«)) n (b.4) 

Now let us consider the most general parity-odd corrections to the stress tensor 
and canonical entropy current, purely based on symmetry grounds. We will then 
explore the restrictions imposed on these corrections from the demand that the en- 
tropy current be divergence free, parity-odd fluid dynamics was recently examined 
in [20] quite generally (including the presence of a global U(l) charge). In fact, we 
can directly take over their analysis by setting the charges to zero; the parity-odd 
terms which are allowed in the neutral fluid are naturally entropy conserving. Here 
we reproduce their result by restricting ourselves to the case of neutral fluids. 
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As in §3 we now work in the Landau frame for this calculation. At first order 
the possible parity-odd corrections for the neutral fluid is given by [20] 15 

nf; arity _odd) = -^^ + x^^, (B.5) 

where 

cr"" = e^u a a v l (B.6) 

is the parity-odd dual of the shear tensor; the corresponding transport coefficient fj 
is referred to as the Hall viscosity Note that we do not consider any parity-even 
corrections to the stress tensor because we know that the shear and bulk viscosities 
have to vanish to ensure dissipation-free behaviour of our fluid. 

The first order equations of motion (A. 4) allow us to eliminate the gradient of 
the entropy or equivalently the temperature in terms of velocity gradients. As a 
consequence of this equation of motion there is only one transverse pseudo-vector 
and one pseudo-scalar at first order which are respectively given by 

U" = e^ p u v u a V a u p , Q = e" up Upd v u p . (B.7) 

Therefore at first order in gradients there are only two possible on shell linearly 
independent vectors; we choose to write these two vectors as 

1 P (B.8) 

V£ = e» up V u u p = -U" - a 

Using these independent vectors the most general entropy current at first order in- 
volving parity-odd terms is: 

J" { s) = J&n + «1 K + «2 V£ (B.9) 

We could have considered adding parity-even corrections to the entropy current in 
(B.9). These however that do not affect the parity-odd analysis; moreover as de- 
scribed in [20] the coefficients of such parity-even corrections would be set to zero 
when we demand that there is no local entropy production. 

Now using the 2 + 1 dimensional version of (3.7) we can directly evaluate the 
divergence of the canonical part of the entropy current in the parity-odd sector, we 
find 



V M = <V n ^ + o P ^ WV J = — *n © O. (B.10) 



The divergence of the remaining two terms are given by 

Vp ( ai V? + a 2 V£) = ^ V£ V pS + ^ V$ VpS + a, V ^ + a 2 

— + — ai ) Qil 

as as 



(B.ll) 



15 



Note that the fl defined in [20] differs from that defined in (B.2) by a minus sign. 
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Thus divergence of the full entropy current is given by 



V /^ J (s) 



' da.2 dT 
Xn + s[— + — oti 
1 \ as as 



QQ (B.12) 



This implies that for zero divergence of the entropy current we must have xn fixed 
in terms of the coefficients appearing in the entropy current 

Since the parity-odd shear tensor a^ u is passive in the entropy analysis, it follows f] 
remains completely unconstrained. 

Therefore we see that the existence of a divergence-free of the entropy current be 
zero fails to constrain either of the two parity-odd first order transport coefficients 
in (B.5). One of these coefficients, the Hall viscosity, does not enter the divergence 
of the entropy current. The other is determined in terms of the possible corrections 
to the canonical form of the entropy current which are arbitrary at this order. 

There are two pieces of frame invariant data in (B.5), that is both C^ v and Gi 
are non-zero when we plug in (B.5) into (4.2) and they are given by 

1 ' B.14 

Comparing this result with the frame invariant data present in the stress tensor 
following from the action (B.l), we conclude that through the action we are only 
able to capture the effect of the vorticity (QP^ U ) term and xn is determined in terms 
of the function s(s) in the action through the following relation 

*> = -(?» + < B15 > 

The Hall viscosity naively seems to be undetermined from the action formalism. 

In passing we note that it may be possible to make modifications to the stress 
tensor (B.3) as has been argued on physical grounds in [15]; the modified stress tensor 
then can captures the effect of the Hall viscosity. We do not see such modifications 
directly from at the level of the action. In any event without any such modifications 
or improvements to the stress tensor, the one which directly follows from the action 
(B.l) captures only transport corresponding to vorticity Q. 
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